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1 Introduction 



In the last few decades, fractional Kolmogorov or equivalently fractional 
Fokker-Planck equations have appeared as an essential tool for the study of dy- 
namics of various complex stochastic processes arising in anomalous diffusion 
in physics [THIGH], finance [8], hydrology [2J, and cell biology [32]. Complexity 
includes phenomena such as the presence of weak or strong correlations, differ- 
ent sub- or super-diffusive modes, and jump effects. For example, experimental 
studies of the motion of macromolecules in a cell membrane show apparent 
subdiffusive motion with several simultaneous diffusive modes (see \'2'2\). 

The present paper identifies a wide class of stochastic differential equa- 
tions (SDEs) whose associated partial differential equations are represented 
by time-fractional order pseudo-differential equations. This connection pro- 
vides stochastic processes whose dynamics correspond to time-fractional order 
pseudo-differential equations. 

Let B t be an m-dimensional Brownian motion defined on a probability 
space (f2, T 1 P) with a complete right-continuous filtration {Tt). A deep con- 
nection between a stochastic process and its associated partial differential 
equation is expressed through the Kolmogorov forward and backward equa- 
tions [JJ. This concept is based on the relationship between two main compo- 
nents: (i) the Cauchy problem 

^ = Au(t, x), u(0, x) = tp(x), t>0,x£ W 1 , (1.1) 
where A is the differential operator 

n d 1 n <9 2 

j=l J i,j=l 

with coefficients bj(x) and (Jij(x) satisfying some regularity conditions; and 
(ii) the associated class of Ito SDEs given by 

dX t = b{X t )dt + (j{X t )dB tl X = x. (1.3) 

The coefficients of SDE (jl.3p are connected with the coefficients of the operator 
A as follows: b(x) = (px(x), . . . , b n (x)) and cri,j(x) is the (i, j)-th entry of the 
product of the n x m matrix cr(x) with its transpose a T (x). 

One mechanism for establishing this relationship is via semigroup theory, 
in which the operator A is recognized as the infinitesimal generator of the semi- 
group T t (-)(x) := ¥,[(■) (X t )\Xo = x] (defined, for instance, on the Banach space 
Cb(R n ) with supnorm), i.e. Aip(x) = lim t _j.o (Tt — I)f(x)/t, if S Dom(^l), the 
domain of A. A unique solution to (|1.1[) is represented by u(t,x) — (T t ip)(x). 

Enlarging the class of SDEs in f| 1 . 3[) to those driven by a Levy process 
leads to a generalization of connection (i)-(ii) where the analogous operator 
on the right-hand side of (jl.lj) has additional terms corresponding to jump 
components of the driving process (see [Tll24| and references therein) . In this 
case, the operator A in (|1.2|) takes the form C(x, D x ) in (|2.14p . 
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A fractional generalization of the Cauchy problem (jl.ll) with A = C(x, T) x ), 
in the sense that the first order time derivative on the left side of equation (jl.lj) 
is replaced by a time-fractional order derivative, has appeared in the framework 
of continuous time random walks (CTRWs) and fractional kinetic theory 
ITT l lT§ p g i En ] . Papers [T l fTS P T p g] establish that time-fractional versions of the 
Cauchy problem are connected with limit processes arising from certain weakly 
convergent sequences or triangular arrays of CTRWs. These limit processes are 
time-changed Levy processes, where the time-change arises as the first hitting 
time of level t (equivalently, the inverse) for a single stable subordinator. 

This paper generalizes the connection (i)-(ii) to time-fractional pseudo- 
differential equations that imply a fractional analogue of Kolmogorov equa- 
tions. First, it establishes the class of SDEs replacing (ll.3[) which is associated 
with the following Cauchy problem: 

u(t, x) = C(x, ~D x )u(t, x), u(0, x) = tp(x), t > 0, x e M", (l.l') 

where D* is the fractional derivative in the sense of Caputo with /3 6 (0, 1) 
(see Section[3]), and C(x, D x ) is the pseudo-differential operator in (|2.14[) . The 
driving processes of the associated class of SDEs are Levy processes com- 
posed with the inverse of a /3-stable subordinator, f3 e (0, 1) ("Theorem 13.51 for 
N = 1). Since such processes are semimartingales, SDEs with respect to them 
are meaningful and have the form in (|3.20[) . A partial result when the driv- 
ing process is either Brownian or Levy stable motion with drift time-changed 
by the inverse of a single stable subordinator is considered in |141I15| without 
specifying the explicit form of the corresponding SDEs. 

More generally, the class of SDEs in the above discussion when the time- 
change process is the inverse of an arbitrary mixture of independent stable 
subordinators gives rise to a Cauchy proble m wit h a fractional derivative with 
distributed orders (see (I2.3[l ) on the left of (jl.l | . namely, 

V^uft^) = C(x,T> x )u(t,x). 

In this case, the time-change process is no longer the inverse of a stable sub- 
ordinator if at least two different indices arise in the mixture. Moreover, SDEs 
corresponding to time-fractional Kolmogorov equations cannot be described 
within the classical Brownian- or Levy-driven SDEs. 

Section [5] of this paper recalls the required auxiliary facts. Section [3] for- 
mulates and proves the main results of the paper, and provides examples, 
including a fractional analogue of the Feynman-Kac formula. Section @] illus- 
trates an alternative technique for establishing the main results in the case of 
Brownian motion. 

2 Preliminaries and auxiliaries 

The fractional integral of order (5 > is 

J"g(t) = jmj ^~ u ) ' 1 9( u ) d ^ 1 > °- (2-i) 
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where r(-) is Euler's gamma function. By convention, J° = /, the identity 
operator, and Jg{t) := J l g(t), the integration operator. The fractional deriva- 
tive of order (3 6 (0, 1) in the sense of Caputo is D*g(t) = J 1- * 9 4zg(t), t > 0. 
By convention, set D* = 4r for f3 = 1. The Laplace transform of D*g is ([5]) 

[DS]( S )= S ^(s)-s /3 -V0+) ! (2.2) 

where g(s) = jC[g](s) — g(t)e~ st dt, the Laplace transform of g. 

Let /i be a finite measure on [0,1]. The fractional derivative with distributed 
orders is the operator (see, e.g., [28} ) 

V fl g(t) = f DPg(t)dfi(p). (2.3) 
Jo 

These operators provide a generalization of fractional order derivatives. The 
mapping (3 i-> D* g(i) is continuous on [0, 1) for a differentiable function g. For 
example, if fi = Ci8 Pl + C 2 Sp 2 , then V^t) = C{Dfcg(t) + C 2 B^g{t). 

If a function ij;(x,£) : K™ x M™ — > C is continuous and satisfies a suitable 
growth condition as |£| — >■ oo, then for u € C^°(]R n ), the operator 

Au(x) = 7 -^[ i>(x,0H0e i(x '°dl;, x€R n , (2.4) 

(/^J JR™ 

where = J R „ u(a;)e — ''"^dx, is meaningful and called a pseudo-differential 
operator with symbol ip(x,£). For properties of pseudo-differential operators, 
see the monographs [TUlfTTll^S] . 

Pseudo-differential operators of interest in this paper are infinitesimal gen- 
erators of strongly continuous semigroups constructed from stochastic pro- 
cesses which are solutions to SDEs driven by a Levy process. Such processes 
are Feller processes, and therefore, they have strongly continuous semigroups 
(see [1]). To this end, we will consider symbols ijj(x,£) which are continuous 
in x e R" and £ £ K n and, for each fixed x, both hermitian and conditionally 
positive definite in £ (see [I] for details). 

A Levy process L t € K n , t > 0, with Lo = 0, is an adapted cadlag 
process with independent stationary increments such that for all e, t > 0, 
lim s _>{P(|.Lt — L s \ > e) = 0. Levy processes are characterized by three pa- 
rameters: a vector b G R n , a nonnegative definite matrix U, and a measure v 
defined on R ra \ {0} such that J min(l, |x| 2 )(i^ < oo, called its Levy measure. 
The Levy-Khintchine formula characterizes a Levy process (as an infinitely 
divisible process) in terms of its characteristic function <?*(£) = e t3> ^, with 

nO = iM-hsU)+ I (e^-l-i^.Ox { \ w \<x){w)) V {dw). (2.5) 

Z JR™\{0} 

The function ^ is called the Levy symbol of L t (see, e.g. [Tll21|.) 

Particularly important for this paper is the class of stable subordinators. 
For /3 G (0,1), a (3-stable subordinator is a one-dimensional strictly increas- 
ing Levy process D t starting at which is self-similar, i.e. {D at ,t > 0} has 
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the same finite-dimensional distributions as {a}/@Dt,t > 0}, and the Laplace 
transform for D\ is given by 



s > 0. (2.6) 



It follows from the general theory of Laplace transforms (see, e.g. [30]) that the 
density /d (t) of D\ is infinitely differentiable on (0, oo), with the following 
asymptotics at zero and infinity [T()ll26] : 

f D ( T ) „ / ( " )2(1 '" > e -d-«(f , t — y 0; (2.7) 
y/2*P(l-P) 

/0 - (T) ~ r(i4 1+< " T400, (2 ' 8) 



Consider an SDE driven by a Levy process 
Jo Jo 



(2.9) 



H(Y s _,w)N(ds,dw) + / / #(Y S -, w)N(ds, dw), 

J\w\<l Jo J\w\>l 

where x S R™, and the continuous mappings & : R" — ► R™, cr : R" -} R nxrn , 
H : R" x R" -> R", and if : R" x R™ -> M™ satisfy the following Lipschitz 
and growth conditions: there exist positive constants C\ and C 2 satisfying 

• \b( yi ) - b(y 2 )\ 2 + \\<r{yi) - <j(y 2 )\\ 2 + I \H(y u w) - H{y 2l w)\ 2 v{dw) 

J\w\<l 

<C 1 \y 1 -y 2 \ 2 ,Vy 1 ,y 2 eR n ; (2.10) 
\H{y,w)\ 2 v{dw) < Ca(l + |y| 2 ), Vy £ R n . (2.11) 

M<i 

Under these conditions, SDE (|2.9p has a unique strong solution Y t (see, [Tll24|). 
If the coefficients b, a, H, and K are bounded, then (T t ip)(x) = E,[ip(Y t )\Yo = 
is a strongly continuous contraction semigroup defined on the Banach space 
Co(R"). Moreover, the pseudo-differential equation associated with the process 
Yt takes the form 

= L{x,Ti x )u{t,xl (2.12) 

at 

where the infinitesimal generator C(x,T> x ) is a pseudo-differential operator 
with the symbol 

*(x,Z)=i(b(xU)-±(E(x)te) (2-13) 
+ / (e < ( G ^'«-l-i(G(x,«;),Ox(|™|<i)H)Kd«'). 

JR"\{0} 
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where G(x,w) — H(x,w) if \w\ < 1, and G(x,w) = K{x,w) if \w\ > 1 
([Dill])- F° r each fixed x G R n , the symbol ^{x, £) is continuous, hermitian, 
and conditionally positive definite [3lfTT]. Using = i(d/dxi, . . . ,d/dx n ), 
the pseudo-differential operator £(i,D x ) has the form 

C(x, D x )v?(x) = i(b(x), T> x )ip{x) - ^{2(x)T> x ,T> x )<p(x) (2.14) 
+ / [<P(% + G(x,w)) -ip(x)- ix{\w\<\){w)(G{x,w),T) x )y{x)]v{dw). 

JR"\{0} 

Here, C(x,T> x ) : C* 2 (R") — >• C (R"), i.e. C 2 (R") c Dom(£(a;,D x )). 



3 Main results and examples 



Let D t be an (J r t )-adapted strictly increasing cadlag process starting at such 
that limj-joo D t = oo a.s. The inverse or the first hitting time process E t of 
D t is defined by E t := inf{r > : D T > t}. The inverse E t is a continuous 
(J- 1) -time- change, i.e. it is a continuous, nondecreasing family of (J-t)-stopping 
times. In fact, for any t, r > 0, we have {£( < t} = {D r _ > i} G J>. Hence, 
by the right-continuity of the filtration (Ft), each random variable E t is an 
(J-t)-stopping time. 

Let Di t t and Z?2,t be independent (J-t)-adapted strictly increasing cadlag 
processes. Then D t — D\ t + -Da,* also possesses the same property, and its 
inverse process E t satisfies ¥(E t < r) = P(D T > t) = 1 - (i^ * F^ ] )(t). 
Here, for fc = 1,2, Fi fc) (i) = P(D fc , T < i) with density f^ k) (if it exists), 
and * denotes convolution of cumulative distribution functions or densities, 
whichever is required. For notational convenience, if a, b > 0, let 

rs=t 



,(2) 



(*) == / n 

's=0 



(i) 



t - s 



dF- 



which, if the density functions exist, can also be written as 
where J is the usual integration operator. 



ds, 



Lemma 3.1 Let D t = c\Di^ -Vc^Di^, where c\,ci are positive constants and 
Dij and Di t o,re independent stable subordinators with respective indices j3\ 
and $2 in (0,1). Then the inverse E t of D t satisfies 



V{E t < r) = 1 - 
and has the density 



F 



(i) 



ClT 



1/01 



* F, 



(2) 



C2T 



1/02 



(*) 



(3.1) 



1 



C2T 



1/02 



(Jfi 



(ih 



ClT 



1/01 



(2) 



C2T 



1/02 



(*) • (3-2) 
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Proof Since D\ T and -D2,t are independent and self-similar processes, 

F(E t < r) = P(£ r > t) = 1 - P(cir 1/fll Di,i + c 2 T X ' h D 2<1 < i)> 

by which (|3.1[) follows. Differentiating (13. 1|) with respect to r yields (13. 2[) . □ 

Lemma 3.2 For any t < oo, the density fE t ( T ) * n P-2P *s bounded, and there 
exist a number (3 £ (0, 1) and positive constants C , k, not depending on t, such 
that 

f Et (r) <Cexp(-fcr^) (3.3) 

for t large enough. 

Routine elementary calculations using (|2.7[) and (|2.8p yield Lemma 13.21 
This lemma can be extended for processes E t which are inverses of stochas- 
tic processes of the form D t = J2k=i c kDk,t, where Dk : t,k = 1,...,N, are 
independent stable subordinators of respective indices j3 k £ (0, 1), and Ck are 
positive constants. 

Theorems 13.11 and 13.21 require the following assumption: {T t ,t > 0} is a 
strongly continuous semigroup defined on a Banach space X with norm || • ||, 
such that the estimate 

\\TM\ < M\\p\\e^ (3.4) 

is valid for some constants M > and cj > 0. For example, if {T t } is a 
semigroup associated with a Feller process, then (|3.4[) is satisfied with uj = 0. 
This assumption implies that any number s with Re(s) > uj belongs to the 
resolvent set p(A) of the infinitesimal generator A of T t and the resolvent 
operator is represented in the form R(s,A) = J e~ st T t dt (see [1]). 

Theorem 3.1 Define the process D t = ^2k=i c k-Dk.t, where D^.t, 1 < k < 
N, are independent stable subordinators with respective indices f3k 6 (0, 1) 
and constants Ck > 0. Let E t be the inverse process to D t . Suppose that 
Tt is a strongly continuous semigroup in a Banach space X , satisfies \3.J$ , 
and has infinitesimal generator A. Then, for each fixed t > 0, the inte- 
gral J Q fE t {T)T T ip dr exists and the vector-function v(t) = J Q fE t (T)T T (p dr, 
where ip € Dom(A), satisfies the abstract Cauchy problem 

N 

C k T>P»v(t) = Av{t), t > 0, w(0) = tp, (3.5) 

k=l 

where C k = cf k , k = 1, . . . , N. 

Proof For simplicity, the proof will be given in the case N = 2. First, define 
the vector-function p(r) = T T (p, where ip £ Dom(^l). In accordance with the 
conditions of the theorem, p(r) satisfies the abstract Cauchy problem 

^-=Ap(t), p(Q) = <p, (3.6) 
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where A is the infinitesimal generator of T T . Now consider the integral J °° /e 4 (r)T T tp dr. 
It follows from Lemma 13.21 and condition (|3.4D that 



/ Et (T)||T T¥ ,||dT<C|M| 



oo, 



(3.7) 



where (3 € (0, 1) and C, k > are constants. Hence, the Bochner integral 
Jo fE t { T )T T <P dr exists for each fixed t > 0. Denote this vector-function by 



v(t) :- 



/ f Et (r)T T (pdT. 
Jo 



(3.8) 



The definition of T t implies v(0) = lim t _>o + fE t (T)T T ip dr = T^ip = ip, ii 
the norm of X. By 



«(*) = 
Since 



i 



o <9r c 2 r 



ClT 



1/01 



(2) 



C2T 1 



(i) > T T tfdr. 



f(2) 



1 1 
6 
1 



Gs) = --(a/ 1 (1) M) 



= - s f[ 1 \as)f?\bs), 



using (12.61) . the Laplace transform of v(t) takes the form 

v(s) = ~ jf° e"^ 2 ^ 2 } 2V„ (3.9) 

= (cf s^+dV 2 " 1 ) [2>](cf a* + c§V a ) 

where — cf fc , fc = 1,2. Due to (|3.4[) . this is well defined for all s such that 
Cis^ 1 + C 2 s^ 2 > w. On the other hand it follows from (j3~6| that 

(s - .A)p(s) = ip, Vs > w. (3.10) 

Let w > be a number such that s > uj q iff Cis^ 1 + C^ 2 > w. Then (|3~^)l 
and (|3.10D together yield 

[Ci^ 1 + C 2 s h - A]v(s) = (Ci^ 1 - 1 + C 2 s^- 1 )p, s > ojo- 

Writing this in the form 



d[s ft t;(s) - s^viO)} + C 2 [s^v(s) 



o/32-l 



v(0)] — Av(s), S > LJ , 



recalling (|2.2p , and applying the inverse Laplace transform to both sides gives 

C l T>^v{t) + C 2 T>^v(t)=Av{t). 

Hence v(t) satisfies the Cauchy problem (13.51) . The case N > 2 can be proved 
using the same method. □ 
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The next theorem provides an extension with D t as the weighted average 
of an arbitrary number of independent stable subordinators. It is easy to verify 
that the process D t — J2k=i c kDk,t given in Theorem 13.11 satisfies 



InEle 



N 

* _1 fe=l 



The function on the right-hand side of (|3.11l) can be expressed as the integral 
— Jq 1 s@dfj,(/3), with fi the finite atomic measure /i = X^fcli c k k ^k- The inte- 
gral fifPdptf) is meaningful for any finite measure p, defined on [0, 1]. Let 
S designate the class of cadlag (J-t)-adapted strictly increasing processes Vt 
whose Laplace transform is given by 

E[e- sVt ] =exp\-t [ s f3 dn(f3)] , s > 0, (3.12) 

where /_i is a finite measure on [0, 1]. By construction, Vb = a.s., and 14 can 
be considered as a weighted mixture of independent stable subordinators. For 
the process Vt £ S corresponding to a finite measure fx, we use the notation 
Vt — D(fi; t) to indicate this correspondence. 

Theorem 3.2 Assume that D(n; t) 6 S where \x is a positive finite measure 
with supp /i C (0,1), and let E t be the inverse process to D(fj,;t). Then the 
vector-function v(t) = / °° fE t (T)T T ip dr, where T t and ip are as in Theorem 
\3.1l exists and satisfies the abstract Cauchy problem 

V M v(t):= f TiPv(t)d(i(P) = Av(t),t>Q,v(0) = (p. (3.13) 



Proof We briefly sketch the proof since the idea is similar to the proof of 
Theorem 13.11 Since supp n C (0,1), the density fD(fi;t)(T~),T > 0, exists and 
has asymptotics (12.71) with some f3 — (3 n £ (0, 1) and (|2.8[) with some (3 — (3\ 6 
(0, 1). This implies the existence of the vector-function v(t). Further, one can 
readily see that v(t) = — J °° J f D{u;r){t)} iT T ip)dr. Now it follows from the 
definition of D(/j,; t) that the Laplace transform of v(t) satisfies 

= JoVW r e -rf^^ m( T T(p )dr = O&ptoW), 3 > O, (3.14) 
s Jo s 

where 77(5) = J Q s /3 c?/i(/3), p is a solution to the abstract Cauchy problem (I3.6p . 
and a) > is a number such that s > Q if 77(5) > u (u> is uniquely defined since 
r)(s) is a strictly increasing function). Combining (|3. 141) and (|3.10[) . 

Ms) - A)v(s) = s > w. (3.15) 

s 

Applying the Laplace transform to (|3.13p yields (|3.15[) . as desired. □ 
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Remark 3.1 

a) Equation (|3.13[) is called a distributed order differential equation (DODE). 

b) If oj = in (13.4[) . that is the semigroup T t satisfies the inequality 
||7t|| < M, then the condition supp /i C (0, 1) in Theorem l3.2l can be replaced 
by supp /i C [0, 1). 

Suppose that L t is an (J 7 t )-adapted Levy process. Let E t be a continuous 
(J r i )-time-change. Consider the following SDE driven by the time-changed 
Levy process L Et ■ 



f b(E s ,X s _)dE s + [ a(E s ,X s ^)dB E 
Jo Jo 



X t = x+ I b(E s ,X s -)dE s + I a(E s ,X s _)dB Es (3.16) 
H(E s ,X a -,w)N(dE a ,dw)+ I I K(E s ,X s _,w)N(dE s ,dw), 

\w\<l Jo J\w\>\ 

with ieR, and continuous maps b(u, y) : R+ x E™ R™, a(u, y) : R + xlM 
R nxm , and G(u,y,w) = X(\ w \<i)(w)H(u,y, w) + X(\ w \>i)(w)K(u, y,w) : R+ x 
R" x R" 4 R" satisfying the Lipschitz and growth conditions ([XTD]) . (ETT|) 
with respect to the variable y S R™, for each fixed u > 0. 

SDE (|3.16p is obtained from an SDE driven by a Levy process upon replac- 
ing its driving process L t by LE t - Since L t is an (J r t )-semimartingale, it follows 
from Corollary 10.12 of 12 that Ls t is an (J^J-semimartingale. Thus, (|3.16[) 
is the integral form of an SDE driven by an (J^J-semimartingale. We use the 
following shorthand form for the differential problem given by SDE Q3.16P : 

dX t = F(E t ,X t ^)QdL Et , X a = x, (3.17) 

where F(u,y) = (b(u,y),a(u,y),G(u,y, •)) indicates the triple of coefficients 
controlling the drift, Brownian, and jump terms, respectively. In the future, 
the integral form in p,16[) will also be expressed using the symbol 0. 

The SDE in ([3T7) is closely related to the following SDE driven by the 
Levy process L t : 

dY T = F(t, F T _) dL T , Y = x. (3.18) 

A duality between SDE (|3.17[) and SDE (13.18[) exists as a special case of a 
general duality result in |13j . but what is required here is given below. 

Theorem 3.3 Let D t be a cadlag (Tt)- adapted strictly increasing process, and 
E t be its inverse. 

1) lfY T satisfies SDE fiTWj} . then X t := Y Et satisfies SDE [3l7\ ). 

2) IfX t satisfies SDE (3~T1\ ), then Y T := X Dt satisfies SDE JO^] . 

Proof Since D t is a strictly increasing (J r t )-adapted process, its inverse E t is 
a continuous (J r t )-time-change. Suppose that Y T satisfies SDE (|3.18l) and let 
X t = Y Et . Then by Proposition 10.21 in [H], 

X t = x+ [ F(s,Y s ^)&dL s =x+ [ F(E s ,Y E(s) _)OdL E .. (3.19) 
Jo Jo 
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X t will satisfy SDE (|3.17p . provided that X s _ = (Y oE) s - can replace Y E r s \_ 
in (|3.19p . The equality Y E / S \_ = (Y o £%_ fails only when s > and E is 
constant on some closed interval [s — e, s) C (0,i] with e > 0. However, the 
integrator L o E on the right-hand side of (|3.19p is constant on this interval. 
Hence, the difference between the two values Y E t s \_ and X s _ = (Y oE) s _ does 
not affect the value of the integral. Consequently, (|3.19p is valid with X s - in 
place of Ye( s )~- Thus, X t satisfies SDE (I3.17p . establishing part 1). Similarly, 
part 2) can be proven using Theorem 3.1 in [15] . instead of Proposition 10.21 
in [12]. □ 

Theorem 3.4 ([T1I24]) If F(u,y) satisfies the Lipschitz and growth conditions 
(12.1 Op and (I2.1ip for each u > 0, then SDE H3.18\) has a unique strong solution 
with cadlag paths. 

Corollary 3.1 If F(u,y) satisfies the Lipschitz and growth conditions (I2.10p 
and (|2.11l) for each u > 0, then SDE \S. 17\ l has a unique strong solution with 
cadlag paths. 



Proof A strong solution to SDE (|3.17p clearly exists by Theorems 13.31 and [ 
To prove the uniqueness, notice that the driving process L Et of SDE (13.171) 
is constant on any interval [D S —,D S ] and so is the solution X t . This implies 
a unique representation X t = Xr) E{t) — Y Et , where Y T is a unique strong 
solution to SDE ([3TT51) . □ 

Theorem 3.5 Let Dk t t, k = l,...,iV, be independent stable subordinators 
of respective indices /3f- G (0, 1). Define D t = Sfc=i c kDk,t, with positive con- 
stants Ck, and let E t be its inverse. Suppose that a stochastic process Y T satisfies 
the SDE (|2.9p driven by a Levy process, where continuous mappings b, cr, H , K 
are bounded and satisfy condition (I2.10p . Let X t — Y Et . Then 

1) X t satisfies the SDE driven by the time-changed Levy process 



X t = x+ f b(X s _)dE s + f a(X s 
Jo Jo 



)dB Es (3.20) 



H{X s -,w)N(dE s ,dw) + [ [ K(X s -,w)N(dE s ,dw); 

J\w\<l Jo J\w\>l 

2) if Y T is independent of E t , then the function u(t,x) = E[ip(X t )\X() = x] 
satisfies the following Cauchy problem 

N 

^2C k T)P k u(t,x)=£(x,T) x )u(t,x), u(Q,x)=<p(x), i>0,i£M", (3.21) 
fc=i 

where ip £ C^K"), Ck = c^ k , k = 1, ...,N, and the pseudo-differential 
operator Clx,!)^ is as in \2.11$ with symbol in &2.13\) . 
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Proof Again, for simplicity, we give the proof in the case N = 2. 

1) Since D t is a linear combination of stable subordinators, which are cadlag 
and strictly increasing, it follows that D t is also cadlag and strictly increasing. 
Hence, it follows from Theorem E3] that X t = Y Et satisfies SDE f^2Ul) . 

2) Consider T Y tp(x) = E[ip{Y T )\Y = x], where Y T is a solution of SDE (|23|) . 
Then T Y is a strongly continuous contraction semigroup in the Banach space 
Co(R") (see pQ) which satisfies p.4[) with to — 0, has infinitesimal generator 
given by the pseudo-differential operator C(x, D x ) with symbol ^(x, £) defined 
in ([2TT5j) . and C$(R n ) C Dom(£(a;, D z )). So the function p Y {r,x) = T Y tp(x) 
with if G Cq(R") satisfies the Cauchy problem 

dpY !£ X) =jC(x,-D x ) P y (t,x), p Y (0,x) = ip(x). (3.22) 

Furthermore, consider p x (t,x) — E[(p(X t )\Xo — x] = E[(p(Y Et )\Yo = x] (recall 
that E = 0). Using the independence of the processes Y T and E t , 

roo poo 

P x (t,x)= / E[ i p(Y T )\E t = T,Y = x}f Et (r)dT= / f Et {T)T Y ^{x)dT. (3.23) 
Jo Jo 

Now, in accordance with Theorem 13.11 p x (t,x) satisfies the Cauchy problem 
(l3~2il . □ 

Theorem 3.6 Assume that Z3(/i;i) € 5, where /i is a positive finite measure 
with supp fiC [0, 1), and let E t be its inverse. Suppose that a stochastic process 
Y T satisfies SDE 12. 9\) . and let X t = Y Et ■ Then 

1) X t satisfies SDE UTM) : 

2) if Y T is independent of E t , then the function u(t,x) = ¥,[ip(X t )\Xo = x] 
satisfies the following Cauchy problem 

V^u(t, x) = C(x, T> x )u(t, x), u(0, x) = ip(x), t > 0, x e M". (3.24) 

Proof The proof of part 1) again follows from Theorem 13.31 Part 2) follows 
from Theorem 13.21 in a manner similar to the proof of part 2) of Theorem 

□ 



Remark 3.2 Theorems 13.51 and 13.61 reveal the class of SDEs which are associ- 
ated with the wide class of DODE pseudo-differential equations. Each SDE in 
this class is driven by a semimartingale which is a time-changed Levy process, 
where the time-change is given by the inverse of a mixture of independent 
stable subordinators. Therefore, these SDEs cannot be represented as classical 
SDEs driven by a Brownian motion or a Levy process. 

Corollary 3.2 Let the coefficients b, a, H, K of the pseudo-differential oper- 
ator C(x,T) x ) defined in {2.1$ with symbol in \2.13]) be continuous, bounded, 
and satisfying condition (|2. 10|) . Suppose ip G Cg(IR n ). Then the Cauchy prob- 
lem 

V^ufax) = C{x,T> x )u{t,x), u(0,x) = ip(x), t>0,xe M", 
has a unique solution such that u(t,x) G Cq(R") for each t > 0. 
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Proof The result follows from representation (|3.23D in conjunction with esti- 
mate (|37f| . □ 

Example 1. Time-changed a-stable Levy process. 

Let L a j be a symmetric n-dimensional a-stable Levy process, which is a 
pure jump process. If p L (t,x) — E[<p(L ai t)\L aiT — x], where <p G C 2 (M. n ), then 
p L (t,x) satisfies in the strong sense the Cauchy problem (\23\) 



dt 



K a (~A) a ^p^(t,x), p^(0,x) = <p(x), t > 0, x G K", (3.25) 



where constant depending on a, and (— A) a ' 2 is a fractional power of 

the Laplace operator. The operator on the right-hand side of (|3.25p can be 
represented as a pseudo-differential operator with the symbol t/>(£) := 
Now suppose that Y t solves the SDE 

dY t = g(Y t -)dL a>t , Y = x, (3.26) 

where g(x) is a Lipschitz-continuous function. Notice that (|3.26p takes the 
form given in (|3.18p with the pure jump process L a ^ as the driving process 
and F(x) = (0, 0, g(x)). In this case, the forward Kolmogorov equation takes 
the form ((23]) 

dpY ^ x) = - Ka (-A) a / 2 {[g(x)] a p Y (t, x)}, t > 0, x G M". (3.27) 

Application of Theorem 13.61 implies that X t — YE t satisfies the SDE 

dX t = g{X t ^)dL a . Et , X Q = x, (3.28) 

where E t is the first hitting time of the process D(/j,; t) described in this the- 
orem. Moreover, if E t is independent of Y t , then the corresponding forward 
Kolmogorov equation becomes 

V^p x (t,x) = -K a (-A) a / 2 {[g(x)] a p x (t,x)}, t > 0, x G E", (3.29) 

where 2? M is the operator defined in (|3.13p . When the SDE in (|3.28l) is driven 
by a nonsymmetric a-stable Levy process, an analogue of (|3.29p holds using 
instead of (|3.27D its analogue appearing in [23] • □ 

Example 2. Fractional analogue of the Feynman-Kac formula. 

Suppose that Y t is a strong solution of SDE (|2.9I) . Let Y G R n be a fixed 
point and q be a nonnegative continuous function. Consider the process Y" t 9 
defined by Y t q = Y t if < t < T g , and Y t q = Y if t > T q , where T q is an 
(J r t)-stopping time satisfying F(T q > t\T t ) = exp(— J* q(Y s )ds) . Then Y q is a 
Feller process with associated semigroup (see [T]) 



(T t V)(y)=E 



exp - / q{Y s )ds <p(y t ) 



Y = y 



(3.30) 



14 



and infinitesimal generator C q (x, D x ) = — q(x) + L(x, D x ), where £(3:, Dj) is 
the pseudo-differential operator defined in (I2.14j) . Let E t be the inverse to a 
/3-stable subordinator independent of Y t . Then it follows from Theorem 13.51 
with N = 1 that the transition probabilities of the process X t — Y Et solve the 
Cauchy problem for the fractional order equation 

x) = \-q{x) + C(x, D x )]u(t, x), u(0, x) = ip(x), (> 0, 1 e 1". 

Consequently, (|3.30p . with X t — YE t replacing Y t , represents a fractional ana- 
logue of the Feynman-Kac formula. □ 



4 Time-changed Ito formula and its application 

This section illustrates a new method of derivation of time-fractional differ- 
ential equations based on the time-changed Ito formula in (13) without using 
the duality principle (Theorem I3.3[) . For simplicity, we consider only the one- 
dimensional case with a Brownian motion as the driving process. Let A* be 
the operator defined as 

A*h(y) = -SL{b(y)h(y)} + ^{a 2 (y)h(y)}. (4.1) 

Theorem 4.1 Let B t be a one- dimensional standard (J- 1 ) -Brownian motion. 
Let D t = *^2k = iCkDk,t, where are positive constants and Dk.t are stable 
subordinators of respective indices fik £ (0, 1). Let E t be the inverse process to 
D t . Suppose that X t is a process defined by the SDE 

dX t = b(X t )dE t + a(X t )dB Et , X Q = x, (4.2) 

where b(y) and o~(y) satisfy the Lipschitz condition (|2.10j) . Suppose also that 
Xn t is independent of E t . Then the transition probability p x (t,y\x) = p x (t,y) 
satisfies in the weak sense the time-fractional differential equation 

N 

^2^U^p x (t,y) = A*p x (t,y), (4.3) 
fe=i 

with initial condition p x (0,y) = 6 x (y), the Dirac delta function with mass on 
x, where A* is the operator in (14. ip . 

Proof For simplicity, the proof is done for N — 2. Let Y t = Xr> t - Then it 
follows that X t = Xr> E . t) = Y^ t , as in the proof of Corollary 13.11 Hence, by 
the independence assumption between Y t and E t , 

/•oo 

p x (t,y)= / p Y (u,y)f Et (u)du (4.4) 
Jo 



in the sense of distributions. 
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Since we are not assuming the duality principle (Theorem I3.3[) . the fact 
that p Y satisfies the classical Kolmogorov equation -^p Y (t,y) = A*p Y (t,y) 
cannot be used here. Instead, we employ the time-changed ltd formula to 
obtain another representation of p x in terms of p Y as follows. Let / £ C£°(R). 
Since X is constant on every interval [D S _,D S ], it follows that X D r s \ = 
Xd s = y s and the time-changed Ito formula in [T3] yields 



f(Xt) ~ f{x) 



f'(Y s )b(Y s )ds 



f'(Y s )a(Y s )dB s 



(4.5) 



E t 



f"(Y s )a 2 (Y s )ds. 



Because / £ C* C °°(R), the process M defined by M u := ft f'(Y s )a(Y s )dB s is 
an (J 7 t )-martingale. Taking expectations in (|4.5p and conditioning on E t which 
has density f Et given in (|3.2I) . we have 



E[f(X t )\X =x}- f(x) 



= / E 



M u + £[f'(Y s )b(Y s ) + ^f"{Y s )a 2 {Y s )Y 



f'(Y s )b(Y s ) + -f"(Y s )a 2 (Y s ) 



Y = x 



E t = u,Y = x 
ds fE t {u)du 



f Et {u)du 



oo pu />oo 



by the assumption that Y t — Xu t is independent of E t . The Fubini theorem 
is allowed since / S C^°(E) and b and a are continuous functions. Using p Y , 
the above can be rewritten as 

E[f(X t )\X =x]-f(x) (4.6) 

\f'(y)Hy) + h"(yW(y)}p Y (s,y)dy ds f Et (u)du 

f(y){ f (JA*p Y {u,y))f Et (u)du\dy, 

where J is the integral operator. On the other hand, reexpressing the left-hand 
side of (I4.6[) in terms of p x yields 



^0 J-oo 



E[f(X t )\X = x] - f(x) = / f{y)p x (t, y)dy - f(x) 



(4.7) 



Since / G C^°(R) is arbitrary and C^°(R) is dense in L 2 (R), comparison of 
(I4.6P and (14.7[) leads to another representation of p x with respect to p Y : 



p x (t,y)-S x (y) 



{JA*p Y (u,y))f Et (u)du 



(4.8) 



in the sense of distributions with p x (0,y) = S x (y). 

Now, we use the two representations (|4.4D and (|4.8p to derive equation (|4.3[) 
with the help of Laplace transforms. The Laplace transform of a function v(t) 
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of the form in (I3.8p . with in (|3.2p . is computed as in (|3.9p . Using this fact 
and taking the Laplace transform of both sides in (I4.4p , we obtain 

(a, y) = (d/- 1 + Cia^-^id^ + C 2 s^,y), s > 0, 

where Cfc = cf fc (fc = 1,2); whereas the Laplace transform of (|4.8p is 

p*(»,y) - -4(y) = (Cx^ 1 - 1 + Cis^- 1 ) JA&VxaP* + C 2 s^,y) 
s 

f!-, sPi- 1 4- f7„aft!-l - — 

A*p Y (C lS ^ +C 2 s l32 ,y) 7 s > 0. 



C lS /3i + C 2 s^ 
Combining these two identities, for s > 0, we have 

C 1 (s^(s,y)-s^- 1 S x (y))+C 2 {s^( S ,y)-s^- 1 6 x (y)) 

= (C lS ft -' + C 2S ft -W(C lS ft +C 2 s^,y) = A^(s,y), 

which coincides with the identity obtained from applying the Laplace trans- 
form to both sides of (|4.3p . □ 

Remark 4-1 If SDE (|4.2p contains an additional term p(X t )dt, then the method 
used in the proof of Theorem l4.1l does not work since the relationship Ye 4 = X t 
does not always follow. Example 5.4 in [T3] yields the following conjecture: if an 
additional term p(X t )dt is included in SDE (|4.2p . where p(y) also satisfies the 
Lipschitz condition, then it is expected that the partial differential equation 
corresponding to (|4.3I) may involve a fractional integral term. 
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